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Abstract 

In this paper, we use the blending functions of Lupa§ type (rational) (p, g)-Bernstein operators 
based on (p, g)-integers for construction of Lupa§ (p, g)-Bezier curves (rational curves) and surfaces 
(rational surfaces) with shape parameters. We study the nature of degree elevation and degree 
reduction for Lupa§ (p, g)-Bezier Bernstein functions. Parametric curves are represented using Lupa§ 
(p, g)-Bernstein basis. 

We introduce affine de Casteljau algorithm for Lupa§ type (p, g)-Bernstein Bezier curves. 
The new curves have some properties similar to g-Bezier curves. Moreover, we construct the 
corresponding tensor product surfaces over the rectangular domain (it, v) £ [0,1] x [0,1] depending 
on four parameters. We also study the de Casteljau algorithm and degree evaluation properties of 
the surfaces for these generalization over the rectangular domain. Furthermore, some fundamental 
properties for Lupa§ type (p, g)-Bernstein Bezier curves and surfaces are discussed. We get g-Bezier 
surfaces for (u,v) £ [0,1] x [0,1] when we set the parameter pi = P 2 = 1. In comparison to g-Bezier 
curves and surfaces based on Lupa§ g-Bernstein polynomials, our generalization gives us more 
flexibility in controlling the shapes of curves and surfaces. 

We also show that the (p, g)-analogue of Lupa§ Bernstein operator sequence Lp n , qn (/, a;) converges 
uniformly to f(x) £ C[0,1] if and only if 0 < q n < p n < 1 such that lim q n = 1, lim p n = 1 and 

n —>-oo n—>oo 

lim p™ = a, lim q™ = b with 0 < a, b < 1. On the other hand, for any p > 0 fixed and p / 1, 

n—>oo n—too 

the sequence L™ (f,x) converges uniformly to f(x) £ C[0,1] if and only if f(x) = ax + b for some 
a, b £ M. 

Keywords and phrases: (p, g)-integers; (p, g)-analogue of Lupa§ Bernstein operators; Limit 
(p, g)-Lupa§ operators; Lupa§ (p, g)-Bezier curves and surfaces; de Casteljau algorithm; Tensor 
product. 

MSC: primary 65D17; secondary J^lAlO, J^1A25, 41A36.: 


1 Introduction and preliminaries 

In 1912, S.N. Bernstein [T[ introduced his famous operators B n : CfO, 1] —>• C[0,1] defined for any 
neN and for any function / £ C[0,1] 

B n(f-’ x )=fl(l ) xk ( 1 ~ x)n ~ kf (jl)' * G Pol¬ 

and named it Bernstein polynomials to prove the Weierstrass theorem [7|. 


( 1 . 1 ) 


In computer aided geometric design (CAGD), basis of Bernstein polynomials plays a significant 
role in order to preserve the shape of the curves or surfaces. The classical Bezier curves [2] constructed 
with Bernstein basis functions are one of the most important curves in CAGD [25] . Apart from this, 
Bernstein polynomials has sevaral applications in approximation theory [7] , geometry and computer 
science due to its fine properties of approximation [T5] . 

In recent years, generalization of the Bezier curve with shape parameters has received continuous 
attention. Several authors were concerned with the problem of changing the shape of curves and 
surfaces, while keeping the control polygon unchanged and thus they generalized the Bezier curves in 

mmm- 


The rapid development of g-calculus [2Sj has led to the discovery of new generalizations of 
Bernstein polynomials involving g-integers m mi m [ 22 ] . 

In 1987, Lupa§ [8] introduced the first g-analogue of Bernstein operator as follows 


( lk]«\ 

n 

) 

k 


fc(fc-l) 


q 2 x k (1 — x) 

Ln,q(f]%) = ^ ' 

fc=0 n {(1 - x) + q^~ 1 x} 

3 =1 

and investigated its approximating and shape-preserving properties. 


n—k 


( 1 . 2 ) 


In 1996, Phillips [2U] proposed another q- variant of the classical Bernstein operator, the so-called 
Phillips g-Bernstein operators and attracted lots of investigations. 


B n , q (f;x) = ^2 
k—0 


n 

k 


n—k —1 

X k n (1 - q s x) 
9 s=0 



X € [0, 1] 


where B n , q : C[0,1] —>• C[0,1] defined for any n G N and any function / G C[0,1]. 


(1.3) 


The q- variants of Bernstein polynomials provide one shape parameter for constructing free-form 
curves and surfaces, Phillips g-Bernstein operator was applied well in this area. In 2003, Orug and 
Phillips [181 used the basis functions of Phillips g-Bernstein operator for construction of g-Bezier 
curves and studied the properties of degree reduction and elevation. 

Recently, Mursaleen et al. Em applied first the concept of (p, g)-calculus in approximation theory 
and introduced (p, g)-analogue of Bernstein operators based on (p, g)-integers. They also introduced 
and studied approximation properties based on (p, g)-integers for (p, g)-analogue of Bernstein-Stancu 
operators, (p, g)-analogue of Bernstein-Kantorovich, (p, g)-analogue of Bernstein-Shurer operators, 
(p, g)-analogue of Bleimann-Butzer-Hahn operators and (p, g)-analogue of Lorentz polynomials on a 
compact disk in [la mi mi mi- 


Let us recall certain notations of (p, g)-calculus . 


For any p > 0 and g > 0, the (p, g) integers [n] p ^ q are defined by 


]p,q =P U +P 71 


-p n 3 g 2 


PQ n ~ 2 + Q nl 


P rl -Q ^^ 
p-q > 

n] q , 

n, 


when p ^ g ^ 1 

when p = 1 
when p = g = 1 


where [n] q denotes the g-integers and n = 0,l,2,---. 
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The formula for (p, g)-binomial expansion is as follows: 


(ax + by)^ q := 2^p“ 

k—0 


-fc)(n 


-fc-1) fc(fe-l) 

q 2 


n 

k 


^n—krk^n—k n k 


{x + y)p, q = (x + y)(px + qy)(p 2 x + q 2 y)---{p n 1 x + q n 1 y), 
(1 - x)% = (1 - x)(p - gx)(p 2 - g 2 x) • • • (p 11 ” 1 - q^x), 


where (p, g)-binomial coefficients are defined by 


n 

k 


Wp,9- 


J p,q W\pd\ n ^]p,g- 

Details on (p, g)-calculus can be found in [51 ITT1 15]. 

The (p, g)-Bernstein Operators introduced by Mursaleen et al. [H] is as follows: 

n—k—1 


Bn lP , q (f,x)= „ ( Li) E p^ L> x k ( p s -q s x) f 


P 


k—0 L J P-9 


fc(fc-i) k 


s=0 


[fc]p,q 


P fc ' n Wp,9 


* G [0,1]. (1.4) 


Note when p = 1, (p, g)-Bernstein Operators given bv 1 1.41 turns out to be g-Bernstein Operators. 
Also, we have (p, g)-analogue of Euler’s identity as: 


n— 1 


(1 - x)l q = J] (p‘ - Q S x) = (1 - x)(p - gx)(p 2 - q 2 x)...(p n ~ 1 - q n ~'x) 


s =0 


= E(-d 


fc (n — k)(n — k — l) fc(fc —1) 


P 


q 


k =0 


n 

k 


Another needed formulae, which can be easily derived from Euler’s identity for | ^ |< 1 is: 


OO k(k — l) i OO 

Y q 2 x _= TT |i + (1 

^{p~q) k [k W 11 1 ^ 


j -1 


r} 


(1.5) 


Again by some simple calculations and using the property of (p, g)-integers, we get (p, g)-analogue 
of Pascal’s relation as follows: 


n 

n—k 

n — 1 

i k 

n — 1 

k 

= q 

P,Q 

k- 1 

+ P 
p,q 

fc 

n 

n—k 

n — 1 

i k 

n — 1 

k 

= p 

p,q 

fc-1 

+ q 

p,q 

fc 


( 1 . 6 ) 

(1.7) 


Motivated by the idea of (p, g)-calculus and its importance in the field of approximation theory 
given by the Mursaleen et al., we construct Lupa§ type (rational) (p, g)-Bezier curves and surfaces 
based on (p, g)-integers which is further generalization of g-Bezier curves and surfaces [H HU jUJ [25]. 


In next section, We present a new analogue, i.e, Lupa§ type (p, g)-analogue of the Bernstein 
functions. 
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2 Construction of Lupas (p, (/)-analogue of the Bern¬ 
stein functions 


We set 


K' n q {t) = 


n 

k 


(n-k)(n-k- 1) fc(fc-l) . j.\n—k 

> 2 Q 2 t K (1 — t) 


.7=1 


( 2 . 1 ) 


and . >bp’q(t) are the (p, ^(-analogue of the Lupa§ g-Bernstein functions |4] of 

degree n on the interval [0,1]. 

Also for substitution u = - f-r where t G [0,1) and u G [0, oo), 


and using Euler’s identity for | ^ |< 1 is: 

oo fc ( fc —!) t, 

E q 2 x*' 

(p - q) k [k\ p , q \ 


n p+ 


k =0 


/c—0 


i-i 


( 2 . 2 ) 


We define 


= 


fc(fc-l) h 
q * u 


(p-q) k [k\ p , q \ n {1 + 


i=o 


j-l 


u} 


(2.3) 


where 6p’|° («) are the (p, (^-analogue of limit Lupa§ g-analogue of the Bernstein functions 19j 
of degree n. 


Theorem 2.1 The Lupa§ (p,q)~ analogue of the Bernstein functions possess the following properties: 
(1.) Non-negativity: b k, ff(t) >0, k = 0,1, t G [0,1]. 

(2.) Partition of unity: 

n 

(*) = !’ *€[0,!]. 


k =0 


(3.) End-point property: 

( 1, if k = 0 

b%(0) = 

\ 0, 

r 

1, if k = n 

^(i) = < 

s 

o 


(4-) (p, q) inverse symmetry: 

b;~ k ’ n (t) = b k A(i-t) = bTk n (t) 

V ’ Q 

for k = 0, 1, . . . , n. 


i—k 

l l 
q ’ p 
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(5.) Reducibility: when p = 1, formula \2. 1\ reduces to the Lupa§ q-Bernstein bases. 


Note: From Euler’s identity for | ^ |< 1 , we have: 

OO k(k—l ) , OO . . j—1 

E f q En , =n M 

t" 0 (p-«) fc [fcU! to VW 

OO 

=>E £’?(*) = 1 ’ * e [0,1)- 

fc=o 

Proof: 

Properties 1,3 and 5 are obvious. Here we only give the proofs of properties 2 and 4. 

Property 2: 

When t = 1, the conclusion is clear; when t ^ 1, we apply the (p,q) analogue of Newtons 
Binomial formula: 

Consider (2) 

Ur 

(n — k) (n — k — 1) fc(fc-l) 


E 

k —0 


n 

k 


P 


-q-^~ r (1 -t) 


i—k 


= E 


A 


fc=0 L J P>9 


(>i-fc)(n-fc-l) k(k-l) „ / t 

P 2 1 2 (1-i) 


1 - 1 


= U(1 - t) + qtj [p 2 (l-t) + q 2 t) .( ^(l - i) + g"- 1 ! 


S— 1 


Hence 


EO')" 1 


k—0 

Property (4) To prove this result, we need following relations: 

fc(2rx —1 —fc) 


[n]p,q = [n]q,p and 


n 


n 

k 


k 

L J 

P,<7 

L J 


(pq)~ 


1 _1 
g ’ p 


(p<?) 


fc(fc-i) 
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Consider 


bZ k ’ n (t) = 


n 

n — k 


(fc)(fc — 1 ) (n-fc)(n-fc-1) 


2 - t n ~ k (1 - t) 


p,q 


n {p 3 Hi -t) + q j H} 

j =1 


n 

n — k 


P 


(fc)(fc-i) jn-kHn-k- 1) . , 

2 q 2 




P 2 q 2 n {^Hr*+ ^jM 1 -HI 

1=1 


n 

fc 


(n-fc)(n-fc-l) (k)(k-l) 

1 2 1 2 j-n—k 


l i_ 

g ’ p 


t u-k ^ _ £)* 


J=1 


lk.n 
3i i 
p ’ g 

in—k,n {,\ 

= & i i (*). 


= &Ti(i-t) 

p ’ g 

in—k,n / 


g ’ p 


The Lupa§ (p, < 7 )-Bernstein blending functions for n = 3 are as follows: 

,0,3 = P 3 (1~H 3 

™ (p(l-i) + «i) (p 2 (l-t) + g 2 t) 

,13 _ (p 2 +M + 9 2 ) Pi(l -H 2 

™ (p(l-t) + 9i) (p 2 (l-t) + g 2 H 

,2 3 _ (p 2 +pg + <? 2 ) gi 2 (! -t) 

P ’ 9 ((p(l-i)+«i) (p*(l-t) + qH)) 


P ’ q (p(l-t)+qt) (j?(l-t) + qH) 




q = 0.7 and p = 0.8 q = 5 and p = 10 


(a) q = 0.7, p = 0.8 


(b) q = 5, p = 10 


Figure 1: Lupa§ cubic Bezier blending functions 


Figure [To] and [Tb] show the Lupa§ (p, g)-Bernstein blending functions of degree 3 for different 
values of p and q . Here we can observe that sum of blending fuctions is always unity. 
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3 Degree elevation and reduction for Lupas (p, q)~ 
Bernstein functions 

Technique of degree elevation has been used to increase the flexibility of a given curve. A degree 
elevation algorithm calculates a new set of control points by choosing a convex combination of the 
old set of control points which retains the old end points. For this purpose, the identities 
and Theorem ED are useful. 

Degree elevation 


q n t 


p n { 1 — t) + q n t 

_ P ^ ~ ^_ b k ' n (t) = ( 

p n (l -t) + q n t p ’ qK> \ 


b k p ’Z(t) = (l - pk+ y +l k] )b k p +^ +1 (t) 


P k [n + l-k}\ jkn+1 


n + 1 


PA 


( t ) 


(3.1) 

(3.2) 


Theorem 3.1 Each Lupa§ (p,q)~ analogue of the corresponding Bernstein function of degree n is a 
linear combination of two Lupa§ (p, q)-analogues of the Bernstein functions of degree n + 1 : 


%%(t) = 


p k [n + 1 — k] 


[n + 1] 


— ' b k ^ +1 (t) + ( 1 - P . . ^ p ' q ) b k+1,n+1 (t) 


p.9 


[n + 1] 


p,q 


(3.3) 


Proof: 


b k ’ n (t) = b k,n (t) (1_ 9 t _|_ < LJl _ 

p ’ q{ ’ p n {l~t) + q n t p n {l-t)+q n t 


p n { 1 — t) + q n t 


q n t 

-\ --- 

p n { 1 — t) + q n t 

Using [3Tl and 13.21 we have 

'p k [n + 1 - fc] 


n 

k 


(re —fc)(n —fc —1) fc(fc-l) 


q "ns— t k (1 - t) 


i—k 


J PA 


n {p j ~ 1 0--t) + q j ~ 1 t} 

1=1 


n 

k 


V 


(re-fc)(re-fc-l) fc(fc-l) 


q-^2— t k (1 -t) 


i—k 


J PA 


bp’q {t) = 


[n + i] 


PA 


PA \ Lfc,n +1 


n{p J ' _1 ( 1 -t) + g j_1 t} 

1=1 


p k+1 [n - k]p, q 


p ’ q w + ( 1 [ n +1] I “P.9 


J PA 




Theorem 3.2 Each Lupa§ (p,q)~ analogue of the Bernstein function of degree n is a linear combina¬ 
tion of two Lupa ( p , q)-analogues of the Bernstein functions of degree n — 1 : 


„n— 1 


(l-t) 


b k ’ n (t) = _-_-_ H_-_i_ 1 _ b k ’ n ~ 1 (t) 

mU pn-qi-q + qn-H p ’ q [> p n - l (l-t) + q n ~H p ’ q U 


(3.4) 
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(3.5) 



pU-kqk-lf 

p™ _1 (l — t) + q n ~ 1 t 


tfc —1,71—1 

u pa 


{t) + 


p n-k-l q k(l _ t ) 

p ra_1 (l — f) + q n ~ 1 t 


b 


k,n— 

PA 


1 


W 


Proof We use the Pascal’s type relations of the (p, g)-Binomial coefficient. 
According to formula 11.71 



fpU-k 

n - 

k - 

- 1' 
-1 

+ q k 

n—1 

k 

\ i) fc(fc-i) , , 

) p 2 q 2 (1 - t) n 

V 



pa 




n 

n -t)+q j ~ i t} 

3 = 1 


or 


b k p ™(t) = 


p n- kq k-l t 


n — 1 

k- 1 




(n-fc)(n-fc-l) (fc —l)(fc —2) 


t 1 *- 1 (1 ~t) 


i—k 


J pa 


p n_1 ( 1 — i) + g ra_1 t 


n—1 

n {P 3 ' -1 ^ -t) +q j ~ 1 t} 

3=1 


pn-k-iqk^ _ t ) 

p n_1 (l — t) + g n_1 t 
p^qK-H 


n — 1 
k 


P 


(n —1 —fc)(n —fc —2) fc(fc-l) 


t k (1-t) 


i—k—1 


n—1 

n {p 3-1 ^ - 0 + 9 3-1 *} 

i=i 


p" 


P £ 1 h fc-l,n-i m , £ 1 U h k,n -1 

1 (1 - i) + W p"- 3 (l - t) + 


w 


or 



(g n - fe 

n - 

k - 

- 1' 
- 1 

+p fe 

n — 1 

k 

\ (n-fcXn-fc-l) fe(fc-l) , , 

) P 2 g 2 (1 — t) n * 

V 



pa 


PA ' 


n 


II (P 3 + 

J=1 




-y-“- ,_£-v- ZL - uk,n 

1 (1 -t) + q n ~ 1 t p ’ q K ' p n ~ l ( 1 - t) + q n ~H p ' q 


(t) 


4 Lupas (p, q )-Bezier curves: 

Let us define the Lupa§ (p, g)-Bezier curves of degree n using the Lupa§ (p, g)-analogues of the Bern¬ 
stein functions as follows: 


P(t;p,g) = £ Pi 


(4.1) 


i—0 


where Pi £ R 3 [i = 0,1,..., n) and p > q > 0. Pi are control points. Joining up adjacent points 
Pi, i = 0,1, 2,..., n to obtain a polygon which is called the control polygon of Lupa§ (p, g)-Bezier curves. 

































4.1 Some basic properties of Lupas (p, q)- Bezier curves. 

Theorem 4.1 From the definition, we can derive some basic properties of Lupa§ (p, q)-Bezier curves: 

1. Lupa§ (p,q)~ Bezier curves have geometric and affine invariance. 

2. Lupa§ ( p,q)-Bezier curves lie inside the convex hull of its control polygon. 

3. The end-point interpolation property: P(0 ;p,q) = Po, P(1 ;p,q) = P n - 

4- ( p , q)-inverse symmetry: the Lupa§ ( p, q)-Bezier curves obtained by reversing the order of the control 
points is the same as the Lupa§ (p,q)-Bezier curves with q replaced by ^ and p replaced by E 
5. Reducibility: when p = 1, formula \4.1\ gives the q-Bezier curves. 

Proof. These properties of Lupa§ ( p , q)- Bezier curves can be easily deduced from corresponding 
properties of the Lupa§ ( p , (^-analogue of the Bernstein functions. Here we only give the proof of 
property 4. 

Let Pt = P n -i, * = 0,1, n, then we have 


P*(f;p,</) = X>i $?(*) 


k—0 


= J2 P *i (!-t) 


k—0 


= P(l-f;-,-). 

p q 


Theorem 4.2 The end-point property of derivative: 

P / (0;p,g) = %f(P 1 -P o ) 

pn i 

P , (1;P>9) = ^=t(P“-P“- 1 ) 

i.e. Lupa§ ( p,q)-Bezier curves are tangent to fore-and-aft edges of its control polygon at end 
points. 

Proof: Let 


P(t;p,?) = X> k &£?(*) = 


E Pk 

k—0 


n 

k 


(n —fc)(n —fe —1) fe(fe-l) , . . , 

r-^- L q “4-^ t k (1 - t) n ~ k 


k =0 


n {p j ~ 1 0 - 0 + q j ~ 1 t} 

3 =1 


w (t- lP , q ) 


or 


p (t;p,q) W(t;p,q) = V(t;p,q) 


then on differentiating both hand side with respect to ‘t’, we have 

P'(t-p,q) W{t;p,q) + P(t;p,q) W '(t;p,q) =V'{t;p,q). 


(4.2) 

(4.3) 
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Let 


then 


Ak(t;p,q) = 


(n — k)(n — k — l) fc(fc-l) 


q 2 t (1 -t) 


n—k 


V(t\p : q) = Y^ Pk Al(t\p,q) 

k =0 

From property 2 of the Lupa§ (p, g)-Bernstein functions, we obtain 

n 

W (t;p,q) = ^2 A k{t-,p,q) 


k =0 


as 


(-4 k(t\P,q) = 


/ _ [n]p,q 


where 


Then 


n — 1 
As — 1 


(n-k)(n-k-l) k(k-l) 


[k\p,q 

[ n ]p,q 
[n ~ fc]p,g 

[n]p,q fc_l, , n -l 


J 

n — 1 

fc 


2 q 2 k r (1 -t) 

(n-fc)Q-fc-l) fc(fc-l) 


n—k 


J P,9 


9 fc - A As - 


q-^~ (n - fc) r (1 - t) 

[ n ]p,4 „n-fe-1 


n—k— 1 


[fcU k-iv—, [n-fc] Pl , 

= A n k ll(t-,p,q) - Dl_ k A n k -\t,p,q) 


p n -«- L (n-k) A n k -\t;p,q) 


/-in _ 

'* “ [fc] M 


[ n ]p,<? „fc~l 




'\P,q m n-k -1 


[n - fc]p,g 




V(0;p,g) = P o p <2 \ W(0;p,<?) = p <2 ' 


hence 


Similarly, we have 


V'(0;p,g) = (C? Pi - P 0 ) P™"'* : 


W (0;p,g) = (C" -£>")p ( 2 \ 


P'(0;p,5) = %f(P!-Po) 

p n 1 


V(l;p, 5 ) = P n q^, W(l;p, g) = 


V'(l;p, 5 ) = (C% P„ - D n x P n _0 2 n ~T~ 2) . 


W'(l;p, g ) = (C£ 
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hence 


P'(l ]P,q) 


^(P n -Pn-1 


) 


Theorem 4.3 Planar Lupa§ (p, q)-Bezier curves are variation diminishing, which means that the 
number of times any straight line crosses the Lupa§ (p,q)~ Bezier curve is no more than the number 
of times it crosses the control polygon. 


Proof. For any polynomial /(f), we denote Z teI g(o,oo) [/(*)] as the number of positive roots of 
/(f) on the interval I. For any vector V = (vo, Vi,v„), we write S~(vo,v±, ...,v n ) to denote the 
number of strict sign changes in the vector V. 

Because the sequence of functions (1, f,... ,t n ) is totally positive on [0,1], then for any sequence 
of real numbers a 0 , a ll ..., a n , Z 0<t< i[a 0 + ait+ +a n t n ] = S~(a 0 + ap+ +a n t n ) < S~(a 0 , ai,..., a n ). 


Let C denote a planar Lupa§ (p, q )-Bezier curve, L any straight line, and let I(C, L) the number 
of times C crosses L. Establish the rectangular coordinate system whose abscissa axis is L. Because 
Lupa§ (p, q)-Bezier curves are geometric invariant, we can denote [xi,yi)(i = 0,1, ...,n) the new 
coordinates of the control points. Let P denote the control polygon and I{P, L) the number of times 
P crosses L. Then we will prove that I(C, L ) < I(P, L). 


We make a parameter transformation. Let u = t £ (0,1), so that u £ (0, oo). Then 


1(6*, L) — Zo <t< i 




L k=0 


= z. 


0<*<1 


E yk 

k =0 


n 

k 


(n — fc)(n — k — 1) k(k- 1) , , s . 

> -^2- L q -V t k (1 - t) n ~ k 


j =i 


= z. 


0 <u<oo 


E yk 

k—0 


n 

k 


P~ 


(n-kHn-k- 1) k(k-l) , 
2 q 2 U 


= Z, 


0<w<oo 


E y k 


k =0 


n { p j ~ i + 

3 = 1 


n 

k 


< S~ 


n 

0 


y o, 


n 

1 


yu 


J p,q 


n 

n 


J P,Q 


= s (yo,yi . ,y n ) 

From l4~3l we can obtain the following two corollaries: 


Corollary 4.4 Convexity-preserving: the planar Lupa§ ( p,q)-Bezier curve is convex, as long as its 
control polygon is convex. 


Corollary 4.5 Monotonicity-preserving: let the control polygon be monotonically increasing (decreas¬ 
ing) in a given direction, then the planar Lupa§ ( p,q)-Bezier curve is also monotonically increasing 
(decreasing). 
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4.2 Degree elevation for Lupas (p, q )-Bezier curves 

Lupa§ (p, q)-Bezier curves have a degree elevation algorithm that is similar to that possessed by the 
classical Bezier curves. Using the technique of degree elevation, we can increase the flexibility of a 
given curve. 


P(t;p,«) = £P k &&(t) 

k—0 


n+1 

p(t,p, q ) = Y,nb k P : n q + \t), 

k—0 

where 


Pk = 11 - 


p k [n + 1 — k) 


P,Q 


k-1 


p k [n + 1 — k] 


p,q 


(4.4) 


L lp,q / \ L lp,q 

The statement above can be derived using the identities (l3.1H and (13.2D . Consider 


P(t-,p,q) = P ^ ^ P (t;p,q) + 

p n { 1 —t)+ q n t 


q n t 


p n (l — t) + q n t 


P (t;p,q) 


We obtain 


P {t;p, q) 



i n + 1 ~ k \ P , q 
i n + 


Pk^ +1 W + 


k=0 


P k+1 [n - k ]p,q 



Now by shifting the limits, we have 


n+1 


P (t;p,q) = Y 


p 


[n + 1 - k] 


[« + 4 


n+1 


“ pSijr'w+E i- 


p k [n + 1 - k\ 


-) P k-i^ +1 W 


fc =0 v 1 J P ’9 7 fc =0 

where P+is defined as the zero vector. Comparing coefficients on both side, we have 

p k [n + 1 - k\ \ f p k [n + 1 - k] 


Pk = I 1 - 


[n + 1] 


p,q 


k-l 


[n + !] 


p,q 


where k = 0,1, 2,..., n+1 and P_i = P n+ i = 0. 


When p = 1, formula 14.41 reduce to the degree evaluation formula of the q-Bezier curves. If we 
let P = (Po,Pi,...,P rt ) T denote the vector of control points of the initial Lupa§ (p, q)-Bezier curve 
of degree n, and P (1 ) = (Pq , P*,..., P* +1 ) the vector of control points of the degree elevated Lupa§ 
(p,q)-Bezier curve of degree n + 1, then we can represent the degree elevation procedure as: 

P (1) = T n+1 P 

where 
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T n +1 — 


1 


[n + 1 \p t q 0 

[n + l]p,q - p[n] Pt g p[n\ p , q 


0 

0 


[n + 1 ] p , g 


0 

0 

0 


■■ [n + lU-r 1 ^,, P n_1 [2k 9 

0 ... [n + l] p , g -p n [l] P ,g 

0 ... 0 


0 

0 


0 

p"[l]p,9 

l n + l]p,g 


For any Z E N, the vector of control points of the degree elevated Lupa§ (p, g)-Bezier curves of 

degree n + l is: P (1 ) = T n+ i T n+2 .T n+ iP. As l —> oo, the control polygon Pd) converges to a 

Lupa§ (p, q )-Bezier curve. 


5 Construction of (p, t/)-analogue of Lupas operators 
and its limit form 

In this section, we present (p, g)-analogue of Lupa§ Bernstein operators and its limit form as follows: 


The linear operators L™ q : C[0,1] C[0,1] 




f p n k Wp,?') 

n 

y Np.9 J 

k 


(n — k)(n — k — l) fc(fc-l) 


q 2 x k (1 — x) 


i—k 


k =0 


II {P j ~ 1 ( 1 ~x) + q j ~ 1 x} 

3 =i 


(5.1) 


is (p, g)-analogue of Lupa§ Bernstein operators. 

Again when p = 1, Lupa§ (p, (^“Bernstein operators turns out to be Lupa§ g-Bernstein operators 
as given in mm- 


It follows directly from the definition that operators L™ (f, t) posses the end point interpolation 
property, that is 


Ll Q (m = f(0), L" (/,!) = /(!) 


for all p > q > 0 and all n = 1,2,. 


Now we show that (p, <?)-analogue of Lupas operator reproduces linear and constant functions. 


By some simple computation, we have 


Some auxiliary results: 


( 2 ) L- q (t,^r 1 ) = ^ ri 

(Q\ T n ( f 2 u \ __ u p , qu / qu 
^p,qK 1 - > u+1 [n] p ,q ' u+l^p+qu 


[ n Lp.q 

[n] P) q 


(n+2)x(n+l) 
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or equivalently for x = 

( 1 )^( 1 ,*) = ! 


i — 1 ™ 2 2 


( 2 ) L p,q(t,x) =X 

(Q\ T n (f 2 \ _ p"" 1 ^ , q 2 a- 2 [n-l] P 
V j p,gv ' > ' [n] Pi9 p(l — x)+qx [ra]p, 9 

Proof: (1) is obvious using (p, g)-analogue of Lupas Bernstein basis function. 

(n — k)(n — k — l) fc(fc-l) 




n 

k 


P 


-q-^- t k (1-ty 


= 1, 


fc=0 


n{p j i (i-t)+^ 

t=l 


( 2 ) 


W^) = E 


fr""‘ 

n 

y [ n ]p,q y 

k _ 


( n-fc)(n-fc-l) fc(fc-l) , , 

) 2 q 2 t (1 — t) 


J P,<7 


fc=o 


= E 

k—0 


= E 

k—0 


= E 

k=l 


Yl{p j 1 0~ t ) + q j M 

j=i 


( P n ~ k 

n 

y [ n ]p,q J 

k 


p 


2 q 2 t K (1 — t) K 


J M 


n—1 

n {yi 1 - 1 ) + gjf} 

i=o 


rji—k 


n — 1 
it- 1 


p 


(n-k)(n-k- 1) fc(fc-l) , 

2 g 2 (1 — i) n ^ 


J P,4 


n—1 

n {p j (! - *) + 9^} 

i=o 


r-zn—k 


n — 1 

jfc-1 


p" 


(n-k)(n-k- 1) fe(fc-l) . 

q 2 u 


n— 1 

n (p j + 

3=0 
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n — 1 P 

k =0 


n—k—1 


n — 1 
k 


P~ 


(n — k — l)(n — k — 2) fc(fc +1) 

2 q 2 u 


fe+i 


n—1 

n {p j +Q iu } 

3=0 


U+ 1 


u+ 1 


„n-l 


E 

fc—o 


E 

fc =0 


n — 1 
k 


(n-k-l)(n-k-2) fc(fc-l) fc 

P 2 9 2 (V) 


n — 1 
k 


n—2 

n {p’p + g- 5 ' (<?«)} 

3=0 

(n fc 1 )(ti jc 2) fc(fc-l) fc 

P 2 <7 2 ('V) 


n—2 

7=0 


u + 1 

or equivalently for a: = 

L p, q (t,x) =x. 

+2 _ p" -1 x , <? 2 a- 2 [n-1] 


(3) To prove, L^,x) = fof + 
Consider 


p,q 

p,q 


^5TT)-E 


P n fc [ fc ]p,<7 


\ 2 

n 

j 

fe 


P" 


(n —fc)(re —fc —1) fc(fc-l) , 

2 q 2 


J 


fc=o 


= E 

k—0 


n—1 

n {p j + q j u} 

3=0 


p 


2n—2k [fc] P ,g 


n — 1 
k- 1 


P" 


(71 — fc)(n — fc — 1) fc(fc-l) , 

- q 2 u 


n—1 

n {p j + 

7=0 


u n—l—k [fc + l]p,q 

u + 1 ^ [«]p,9 


n — 1 
k 


(fr 


n—2 

n {p 7 '+9 j (fw)} 

7=0 
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^-r) = 


U [n - 1 ]p, q ^ ( p k + q[k]p,c 


u+1 u+1 [n\ P}q — 


E 


[«-!] 


P.9 


n — 1 
k 


(n — fc — l)(n— fc — 2) fc ( fc !) / qu \k 

V 2 q 2 (y) 


J p.g 


n—2 

n +9 j (f«)} 

j=o 


u p 


n—1 


Z( + 1 [n] Pi q u+1 [n] 


gu [n-l] Pl ,^4 


n — 2 

k-1 


(n — fc — l)(ra — fc —2) fc(fc-l) 

P 2 9 2 (V) 


J p,<j 


n—2 

n {p 7 '+9 j (fw)} 

J=0 


1/, p 


gu [n - 1 } Pt q — 


<?« n—2 


U+1 [n\p t q U+1 [n}p, q 1 + 


qu 


n — 2 

k 


(n — k — 2)(n—k — 3) fe(fc-l) 

p 2 q 2 ( 


J P.9 


p fc=o 


n—3 

p n-2 J] {pj +gi(^u)} 

J '=0 


u p 


qu 


qu \ [n - l] p , 


P,<7 


u + 1 [n] Pl g u+1 \p + guy [n]p, g 

Theorem 5.1 Let 0 < q n < p n < 1 such that lim p n = 1, lim q n = 1 and lim p™ = a lim g™ = b 

n—too n—too n—too n—t oo 

with 0 < a, b < 1. Then for each f G C[0,1], L” ? (/; x) converges uniformly to f on C[0,1]. 

Proof: Proof is obvious using the following Korovkin’s theorem. 


Let (T„) be a sequence of positive linear operators from C[a,b] into C[a,b]. Then 
lim„ || T n {f,x) - f(x)\\ C [a,b] = 0, for all / G C[a,b] if and only if lim„ ||T„(/ i ,a;) - fi(x)\\ C [ a ,b] = 0, for 
* = 0,1, 2, where f 0 (t) = 1, /i(f) = t and f 2 it 2 ) = t 2 . 


Remark: 

For q G (0,1) and p G (g, 1] it is obvious that lim [n] P n = 0 or ——. In order to reach to 

n—>-oo ’ P 9 

convergence results of the operator L™ (f; x), we take a sequence q n G (0,1) and p n G (g„, 1] such that 
lim p n = 1, lim q n = 1 and lim p" = a, lim g™ = b with 0 < a, b < 1. So we get lim [n ] Pn ,q„ = oo. 

n—>oo n—>oo n—> oo n—t oo n—»oo ’ 


(p, g)-analogue of the limit Lupa§ Bernstein operators 

First, let 0 < q < p < 1 and lim - , , = 1 — ( a ) fe 

We define linear operators L™ q defined on C[ 0,1] as 

! o° 

E /(l-(f) fc )^, if x G [0,1) 

fc=0 (5.2) 

/(1), if x = 1, 

which is (p, g)-analogue of the limit Lupa§ Bernstein operators. 

For more details on Limit Lupa§ g-analogue of Bernstein operators, one can refer mm- 


Note that the function L™ {f\x) is well-defined on [0,1] whenever f(x) is bounded on [0,1]. 


Theorem 5.2 Let f G C[0,1], g(x) = /(1 — x). Then for any p > q > 0, 


L£,(/;i) = L? ifol-i), 


for t G [0,1] 
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Proof: 

The proof of above theorem follow easily along the lines of m and using the following relations, 


[n\ P , q = (pq) n 1 [n]i i i, 


p k [n — k\ p ,q P k "Ml,4 


= 1 - 


q ’ p 


and 


n 


n 

k 

p,q 

_ k _ 


(pqY 


fc(2n —1 —fc) 


1 1 
q > p 


(pq) 


fc(fc-i) 


Remark: For p = 1, this equality coincides with formula (10), Theorem 3 in 


Corollary 1: Let 1 be fixed, / £ C[0,1], and g(x) = /(I —a:). Then, for x £ [0, 1], L” g (/;t) 
converges uniformly to Lff q (f: t) for any p > 0 being fixed where 


p>q 




if0<g<p<l 

L fi (/;*)> if P > 9 > 1 

g 5 p 

An explicit form of the limit function for 0 < q < p < 1 is given by 15.2 


(5.3) 


Theorem 5.3 Let p > g > 0, p / 1 be fixed and f £ C[0,1]. Then L^ q (f: t) = f(t) for all t £ [0,1] 
if and only if f[x) = ax + b for some a , b £ R. 

Corrollary 2: Operators (/; x) reproduces linear functions, that is L™ (at + b;x) = ax + b 
for all p > q > 0 and all n = 1,2,.... where a, b £ M. 

Note: (p, g)-analogue of Lupas Bernstein operator have an advantage of generating positive 
linear operators for all p > q > 0, whereas (p, g)-analogue of Bernstein polynomials introduced by 
Mursaleen et al. generate positive linear operators only if 0 < g < p < 1. 


6 (p, q )-de Casteljau algorithm for Lupas Bezier 

curves : 

Lupa§ (p, g)-Bezier curves of degree n can be written as two kinds of linear combination of two Lupa§ 
(p, g)-Bezier curves of degree n — 1, and we can get the two selectable algorithms to evaluate Lupa§ 
(p, g)-Bezier curves. The algorithms can be expressed as: 

Algorithm 1. 


P?(i;p,,)sP?sP, i = 0,1,2. 


« = CM + y-Ci)‘4U rr‘(‘;p,9) 

r = l, ...,n, * = 0,1,2. ,n — r., 


( 6 . 1 ) 
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or 


P?(i;p, ? )sP?sPl i = 0,1,2. 


Pf«;M) = p f«( f ;p.9) + p r‘ 

r = l,...,n, * = 0,1,2. ,n — r., 


( t',P, q ) 


( 6 . 2 ) 


Then 


n—1 

P(i;p,9) = ^2Pl(t;P,q) = • i(t;p,q) b^ n ~ r (t) = ... = P£ (t-,p,q) (6.3) 

i—0 

It is clear that the results can be obtained from Theorem ( 13 . 21 ) . When p = 1 , formula ( 16 . 11 ) and 
( 16 . 21 ) recover the de Casteljau algorithms of classical g-Bezier curves. Let P° = (Pq, Pi, ..., P n ) T , 
P r = (Pq,Pi, ...., Pn- r ) T , then de Casteljau algorithm can be expressed as: 


Algorithm 2. 


P r (f;p,g) = M r (t;p,q)....M 2 (t;p,q)M 1 (t;p,q) P° 
where M r (t',p, q) is a (n — r + 1) x (n — r + 2) matrix and 


(6.4) 


M r (t;p, q) = 


or 


M r (t;p,q ) = 


p n ~ r (i-«) 

0 


0 

0 


p" r (i-0 

p n-r(l_() + ,»-r( 

0 


0 

0 


r t 

p"- r '(i-0+9"- , 'i 

p"~ r (i-«) 

p r*-r.(l_ t ) +9 n-r t 


0 


p” 

p n-r(l_ t ) +9 n-r t 

P "~ r ~ 1 g(l~*) 

(l-t) + 9 "- r t 


0 


q n r l 

p n- r( 1 _ t ) + g n-r t 

p n-r( 1 _ t ) + g n-r t 


p” r 1 qt 
p Tl ~ r {l- t )+q rl - r t 

p 

P n - r ( 1 -t)+9" — 


0 

0 


p n-r(l_ t ) +g n-r t 

P n ~ r (l-t) 

pn — r( 1 _ t ) +g n-r t 


0 

0 


pq n r t 
p n ~ r (l-t)+q n -rt 

q n r (i -t) 

p n ~ r (l- t )+qn-r t 


0 

0 

0 

q n ~ r t 

p n - r {i-t)+ q n ~ r t_ 

0 

0 

0 

q n ~ r t 

P n ~ r (l-t)+ q n - r t_ 


7 Tensor product Lupas (p, g)-Bezier surfaces on 

[0,1] x [0,1] 

We define a two-parameter family P(u, v) of tensor product surfaces of degree to x ro as follow: 


P{u,v) = b P™qM) b J p^ q2 ( v ), (u,v) G [0,1] x [0,1], (7.1) 

i =0 j =0 

where P ij € R 3 (i = 0,1,..., to, j = 0,1,..., n) and two real numbers pi > qi > 0, p 2 > <?2 > 0, 
6 p ’i™i(u), b{;™ q2 (v) are Lupa§ (p, g)-analogue of Bernstein functions respectively with the parameter 
Pi, q\ and P 2 ,q 2 - We call the parameter surface tensor product Lupa§ (p, g)-Bezier surface with degree 
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to x n. We refer to the Pi.j as the control points. By joining up adjacent points in the same row or 
column to obtain a net which is called the control net of tensor product Lupa§ ( p , g)-Bezier surface. 


7.1 Properties 

1. Geometric invariance and affine invariance property: Since 

m n 

EEfcw <&>) = i, (7-2) 

2=0 J =0 

P (u,v) is an affine combination of its control points. 

2. Convex hull property: P(it,u) is a convex combination of P* j and lies in the convex hull 
of its control net. 

3. Isoparametric curves property: The isoparametric curves v = v* and u = u* of a tensor 
product Lupa§ (p, g)-Bezier surface are respectively the Lupa§ (p, q)-Bezier curves of degree to and 
degree n, namely, 

m / n \ 

PKO = E fc(“) - [0,1]; 

2=0 ' j =0 ' 

n / m \ 

P(u*,v)=^(^P i j^ qi (u*))b^ q2 (v) , v e [0,1] 

j '=0 ' i =0 ' 

The boundary curves of P(u, v) are evaluated by P(u, 0), P (u, 1), P(0, v) and P(l, v). 

4. Corner point interpolation property: The corner control net coincide with the four 
corners of the surface. Namely, P(0, 0) = Po,o, P(0,1) = Po,n, P(l, 0) = P m ,o, P(l, 1) = P m>n , 

5. Reducibility: When pi = p-> = I. formula m reduces to a tensor product g-Be'zier patch. 


7.2 Degree elevation and (p, q)-de Casteljau algorithm for 
Lupas Bezier surface 

Let P(u,u) be a tensor product Lupa§ (p, q)-Bezier surface of degree to x n. As an example, let us 
take obtaining the same surface as a surface of degree (m + 1) x (n + 1). Hence we need to find new 
control points P* • such that 


ra+1 n+1 


p(«.«) = EE p yfe(“) Cw = EE p * 


n+l / 


2=0 j =0 

Let Oi = l- pi +- = 1- P = [n+1 ^ 

1 [m+1] ’ 


2=0 j— 0 


[«+iL 


Then 

P tj = Oii fa Pi_i,j_i + at (1 - Pj) Pi-i,j + (1 - ai) Pj P i,j—i + (1 - <+) (1 - Pj) P itj 

which can be written in matrix form as 

pj KH-j] P2 ,„ 


1 _ pj pj [m+1—4 P1 ,„ 

[m+l]_ ~ [ra+ll ” 

L Jpi.qi L J Pl,qi 


Pj— l,j— 1 Pi—l,j 
Pi.i—1 Pi, .7 


1 - 


[«+!]„ 


P2 1^+!—J'lpg.q 
N+dp,.,, 


(7.3) 


(7.4) 
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The de Casteljau algorithms are also easily extended to evaluate points on a Lupa§ (p, q )-Bezier 
surface. Given the control net Pj j £ R 3 , * = 0,1,to, j = 0,1, n. 


p !j°( w - *■’) = P°j° = p iJ * = 0,1,2., to; j = 0,1,2...n. 


P'f(«,«)= 


p? G 1 -") 


h 


pr-l,r-l pr-l,r-l 

i,j 

j-»r—l,r—1 -p»r—l,r—1 

*i+l,j *i+l,j+l 


_p? r (i-“)+«r r “ p r r { i - u )+<ir 
, r=l,...,fc, fc = min(m, n) * = 0,1,2 ., to — r; j = 0,1, ....n — r 


P?~ r ( l-») 


<?" r v 

P2~’'( 1 - v )+ ( l 2 


(7.5) 


or 


P°f(u,v) = P°f = P itj i = 0,1,2 ., to; j = 0,1,2... 


n. 


p-’/M = 


pT ' ’’gi( 1 ~ M ) 


Pi 


<h u 


pr-l,r-l pr-l,r-l 

i,j 2,^+1 

T-jf—1 ,7*—1 Tjr—1,7*—1 

*i+l,j+l 


_p? r (i-«)+ 9 r r “ p™ r ( i - u )+9r r 
, r=l, ...,fe, fc = min(m,n) * = 0,1,2 ., m — r; j = 0,1, ....n — r 


pr j - r qj(. i-G 

Pj - ’'! 1 -' u )+'?2 _r ' i 


.p? r (l—ul+qj r i). 


(7.6) 

When to = n, one can directly use the algorithms above to get a point on the surface. When to **, 
to get a point on the surface after k applications of formula JL3 or (17.61) . we perform formula (16.41) 
for the intermediate point P t . 

Note: We get Lupa§ q-Bezier curves and surfaces for (u,v) £ [0,1] x [0,1] when we set the 
parameter pi = P 2 = 1 as proved in [4j . 
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8 Shape control of (p, q )-Bezier curves and surfaces 

We have constructed Lupa§ type (p, q)-Bernstein functions which holds the end point interpolation 
property as shown in following figures of curves and surfaces. It can be also observed that the curves 
(surfaces) generated is contained within the convex hull of the control polygon (control net) for 
different values of p and q. 

Parameter p and q has been used to control the shape of curves and surfaces: ifO<g<p<l, 
as p and q decreases, the curves (surfaces) move close to the control polygon (control net), as p and 
q increases, the curves (surfaces) move away from the control polygon (control net). 

If 1 < q < p, the effects of p and q are opposite, as p and q decreases, the curves (surfaces) move 
away from the control polygon (control net), as p and q increases, the curves (surfaces) move close to 
the control polygon (control net) which can be seen in the following figures. 




(a) 0 < q < p < 1 (b) 1 < q < p 



Figure 2: The effect of the shape of Lupa§ (p, q )-Bezier curves 
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(c) qi = Q 2 = 3, Pl = P 2 = 5 


(d) qi = q 2 = 20, p 1 = p 2 = 200 


Figure 3: The effect of the shape of Lupa§ (p, g)-Bezier surfaces 


Acknowledgement: We are very grateful to the anonymous referees for the inspiring comments 
and the valuable suggestions which improved our paper considerably. 
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